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Abstract 

In this paper, the problem of finding an axisymmetric stationary space- 
time from a specified set of multipole moments, is studied. The condition 
on the multipole moments, for existence of a solution, is formulated as a 
convergence condition on a power series formed from the multipole mo- 
ments. The methods in this paper can also be used to give approximate 
solutions to any order as well as estimates on each term of the resulting 
power series. 

1 Introduction 

The concept of multipole moments for curved spacetimes has been around for 
more than 35 years, but there are still some unanswered questions. A geomet- 
rical definition of multipole moments for static asymptotically flat spacetimes 
was given by Geroch [5], and was generalised to the stationary case by Hansen 
[TU] . If one also incorporates Beig's pQ generalisation of centre of mass, one get 
a completely coordinate independent definition of the multipole moments. 

One of the remaining questions to answer is whether, for any given set of 
multipole moments, there is a solution to Einstein's vacuum field equations 
with those multipole moments. There are solution generating techniques, for 
instance, for the axisymmetric case, the FHP-method [B] can be used. However 
whether these methods always converge is still an open question. In the general 
static case recent work has been done by Friedrich [7]. However, he establishes 
results for entities related to the multipole moments, but not for the multipole 
moments themselves. In the static axisymmetric case the complete convergence 
conditions were established in [2j- The convergence problem for the stationary 
axisymmetric case is treated in this paper. 

The tensorial recursion of Geroch and Hansen can be simplified to a power- 
series expansion of a function formed from the metric components. This was 
shown in [Sj for the stationary axisymmetric case and in [1] for the general 
stationary case. It is easy to see [1] that a necessary condition on the multipole 
moments, for existence of a solution, is that this power-series expression has 
positive radius of convergence. In this paper, this condition is also proven to be 
sufficient for the axisymmetric case. 
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2 Field equations and potentials 



In this paper we restrict our selves to stationary vacuum spacetimes. In this 
section and the next we will treat the general stationary case but will later 
restrict our selves to the axisymmetric case. 

Consider a stationary spacetime [M,g a b) with time-like Killing vector field 
£ a . Let V be the 3-manifold of trajectories of £ a , i.e., M with time factored 
out. Wc let A = — £ a £ a be the norm of £ a . Furthermore, we define the twist 
uj a through oj a — e a bcd£, b ^7 c £, d - Due to the fact that we only consider vacuum 
solutions, we can define a twist potential uj via V Q o; = uj a , since V[f,w a ] = in 
vacuum. 

The metric g a b (with signature (—,+,+,+)) on the spacetime induces the 
positive definite metric 

h a b = ^9ab + £a£fc 

on V. 

In the stationary vacuum case, Einstein's field equations induce field equa- 
tions on V (rewritten from [9]), 
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Rab = —t((D a \)D b \ + (D a u)D b u) (la) 



D a D a X = -((D a X)D a X - {D a cj)D a uj) (lb) 
A 

D a D a u = j(D a X)D a uj, (lc) 

where R a b and D a is the Ricci tensor and the covariant derivative operator with 
respect to the metric h a b- We introduce the potential 

1_A-^ 

v 1 + A + iw K ' 

With this potential, the field equations reduce to 
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Rab 



— {D {a <f>)D b) <l> (3a) 

D a D a <j> = -3^-(D a <f>)D a <f>. (3b) 

<P<P- 1 

The equation (|3bj) is a useful form of the Ernst equation [5] . The non-linearity 
of this equation makes it quite difficult to solve in general. In this paper however 
we will prove existence of solutions when the multipolc moments are specified. 



3 Definition of multipole moments 

In this section, the definition of multipole moments given by Hansen in [10j . 
is quoted. For the definition one demands that V is asymptotically flat in the 
following sense: There exists a 3-manifold V and a conformal factor satisfying 

(i) V = V U A, where A is a single point 

1 Defined in a region who's spatial part is R 3 minus a ball. 
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(ii) h ab = fl 2 h ab is a smooth metric on V 

(iii) At A, tt = 0, D a n = 0, D a D b n = 2h ab , 

where D a is the derivative operator associated with h ab . Note that we only de- 
mand asymptotic flatness on V. We do not demand that the spacetime "asymp- 
totically tends to Minkowski spacetime". Thus in this sense the Taub-NUT 
solutions are asymptotically flat. 

On M, and/or V one defines the scalar potential^ <j> as in Q . The multipole 
moments of M are then defined on V as certain derivatives of the scalar potential 
4> = (f>/yQ at A. More explicitly, following [ID] , let R ab denote the Ricci tensor 
of V, and let P = <p. Define the sequence P, P ai , P ai a 2 , • ■ ■ of tensors recursively: 

P — C\r> P (n-l)(2n-3) p p i f & \ 

where C[ • ] stands for taking the totally symmetric and trace- free part. The 
multipole moments of M are then defined as the tensors P ai ... a „ at A. 



4 Axisymmetry 

To simplify the situation we will henceforth only consider axisymmetric space- 
times. In this situation there is a canonical way (see for instance |12j ) of finding 
the Weyl-coordinates (t : ip, R, Z), where {-§i) a = £ a and {-§^) a are the time-like 
and space-like Killing vector fields respectively. 
In these coordinates the metric becomes 

ds 2 = -X(dt - Wdip) 2 + X-^Ffdtp 2 + e 2f3 (dR 2 + dZ 2 )), (5) 

where the functions A, [3, and W only depend on R and Z. We also observe 
that A is the norm of the time- like Killing vector field £ a . 
The metric on V becomes 

h ab ~ R 2 dcj) 2 + e 2f3 (dR 2 +dZ 2 ). (6) 

In these coordinates, equation (|3b| reads 



<t>RR + ^ZZ + ^ = (<f R + <P%) • (7) 

Introduce new coordinates p = B ?+ Z i , z = R 2+ Z 2 and rescale the potential 

f=-JL=. (8) 

Vp + z 

To be consistent with the notations in [3], we introduce different conformal 
factors. An initial conformal transformation is formed by fl = (p 2 + z 2 )e _/3 . We 
then get the rescaled metric of the desired form 

h ab = {l 2 h ab - p 2 e- 2/3 d(p 2 + dp 2 + d~z 2 . (9) 



2 This definition differs from Hansen but is equivalent. See theorem 4 in [11] . or for the 
axisymmetric case [3]. 
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The additional freedom in the choice of conformal factor can be described by 
letting fl = e K Q, where k(A) = 0. The vector k'(0) then corresponds to a trans- 
lation. The higher order terms of k leaves the multipole moments unchanged 
and can therefore be freely specified. This freedom was used to simplify the 
recursive definition of the multipole moments in [3j. The metric corresponding 

to ft is h a b = Q 2 h a b = e 2K h a b The potentials are rescaled as 4> — <f>/\/fo, an d 
4> = <j>/y/ti. We note that (f> = e^/ 2 ~ K / 2 ; . = e^ 2 ""/ 2 /- 

The axisymmetric condition also simplifies the description of the multipole 
moments. We find [5] that at A, P ai a 2 ...a„ is proportional to C[z ai z a2 . . . z Q J, 
i.e., we can define scalars m n via 

Pa (A) = m n C[ }, (io) 

where z a = (dz) a is the direction along the symmetry axis. Thus we may 
represent the multipole moments with m n . 

Remark. As well as different sign conventions some authors use the convention 
M n = i? 1 . . . 2 a »P 0l ...on(A). This implies that M n = ±g4m„. 



5 Axisymmetric field equations 

For further studies we will need to reformulate the field equations. The new 
formulation is chosen such that the potential on a coordinate axis is the leading 
order part of the potential. This is useful because as we see later one can 
reconstruct the leading order part of the potential from the multipole moments. 
The definition of the leading order part follows at the end of this section. 

Henceforth we assume that / is real-analytic in terms of (x = pcostp, y = 
psmcp, z) and independent of (p. Thus we only search for solutions of this form. 
These properties implies that / is even in p, where the map p — > — p is treated 

as If — > Ifi + 7T. 

Now define a complex variable w — ^fe^. Note that / will not be analytic 
in w but will have the form 



oo oo 



/H =EE a n , m w n w m . (11) 



n— m— 



The p-symmetry implies the symmetry f(w) — f(w), thus a rhm = a m ,, ; . 
Expressed in w and w the equation (|7|) reads 



f w . = l^LZJl _ V(f + ^f w )(f + ^U) 
2(w — w) 1 — Awwff 

Now we can define a function g to be the function of two independent real 
non-negative variables £, C defined by 

oo oo 
n— m— 

where a„ iTO is the same as in (fTTj) . If ([TTj) converges for \w\ < r then (fT3|) 
converges for < C < ?*o and vice versa. This gives a bijection, mapping / 
to g. Note that the coefficients of / is a m ^ n = a n)m , thus / is mapped to g. 
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The symmetry f(w) — f(w) implies the symmetry <?(£, C) — diCiQ- Hence the 
equation (JT^J) is equivalent to 



g ' c = W^T) ^C9~9 ■ (14) 

Definition 1. Suppose that f is a real analytic function of (x — p cos ip,y = 
psimp,z) independent of (p. We then define the leading order part of f to be 
/i(r) = g(r, 0), where g(£,C) is the function formed from f as above. 

Note that this is equivalent with the definition in [3J . There we treated / as 
a function of (z,p), complexified p via the power-series expression and defined 
h(r) = f(r,-ir). 



6 Specified multipole moments 

In this section we study the relation between the multipole moments and the 
leading order part of the potential /. How to compute the multipole moments 
from the leading order part of the potential and the metric component /3, was 
studied in [3J . Therefore we cite the main theorem from [3J . 

Theorem 2. Suppose that M is a stationary axisymmetric asymptotically flat 
spacetime, with V being a conformal rescaling of the manifold of time-like Killing- 
trajectories, where the metric h a b on V has the form @ . Let <fi be the confor- 
mally reseated potential. Furthermore, assume that <fi and (3 are analytic in a 
neighbourhood of A £ V . Choose n such that, for some constant C, 

KL (r) =-ln(l-r / , dr' - rC) +f3 L (r), (15) 

Jo r 

thereby inducing the metric h a b = e 2K h a b. Define r(p) implicitly by p(r) — 
Put y(p) = 4> L {r{p)) = e- KL{r{ -P^/ 2 4> L {r(p)). Then the multipole 



re 

moments mo, mi, ... of M are given by m n = g-n(Q) 

Note that if we demand / and /3 to be analytic in (x, y, z) and independent 
of ip, the analyticity conditions of the theorem is fulfilled. 

To be able to find a relation between the multipole moments and fu alone, 
we need to eliminate 0l from theorem [H With the choice of Weyl-coordinates 
this can be done as follows. As in [3] we let 

An explicit computation of the Ricci tensor for the metric (0) implies 

VVRab = —, ( 17 

w(w — w) ow 

From (|5a|) we get 
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Combining this with ([5]) gives 

d/3 = w-w _ 2w—)(f I 2w d ^) (19) 

dw (1 — 4wwff) 2 dw dw 

The map / to g can also be done for j3 to get a function &(£, Q. For this function 
we get 

f = IT J i^ + 2 #* + 2 # < 2 »» 

This equation can then be used to reconstruct b and hence (3 from g. 
The observation /3z,(r) = b(r, 0) and /i(r) = <?(r, 0) implies 

^W = (§), e =r/ L (r)+2r/i(r) 2 . (21) 

This is easy to integrate to get (3l and substitute this in theorem^ The constant 
of integration is chosen such that /3l(0) = 0. This must hold due to lemma 3 
in [3], where it was proven that /3 contains a factor p 2 , thus b contains a factor 
(£-C) 2 , i.e., /3 L (r) = C(r 2 ). 

Now we study the problem of finding the potential / given the multipolc 
moments m n £ C. Form 

oo 

71 = 

We will assume that this power series converges in a neighbourhood of the origin. 

To be able to express Jl in terms of y(p) we need the following relation 
between p and r. 



Lemma 3. // y is given by (|22[) and p, k are given implicitly by the relations 
in theorem® then 

= - P -+p f 4 f 2p"\y(p") + 2pW(p")\ 2 dp"dp' + l + K '(0)p, (23) 



r 



o P' 2 



where the variables are real. 

Proof. Theorem [5] gives the relation 



y{p) = e 2 f L ( r (p)) = ^-ff L (r(p)). (24) 



Combining this with (|2ip . we have 



A = ^|y + 2P2/ P | 2 . (25) 



Together the definitions of p and gives 



e/fe-Kt = 1 _ r / d r ' _ r C, 



P ./o r' 2 

dividing by r and differentiating with respect to r gives 

e 2& = rV/p 2 (26) 
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and a further differentiation with respect to p gives the useful relation 

dp p r v p r 

Using this in ([25]) , we get 



2p\y + 2py p \ 2 = -^-(e~ 2 ^). (28 ) 

So that 

e -2/9 t (r) = _ I" 2p \ y + 2pyp \2 df) + Ci (2Q) 
JO 

As proven in [3J, /3l(0) = 0, hence the limit r — > gives Ci = 1. Substituting 
([2"6| and dividing by p 2 gives 

^ - = -7/0 2pW) + 2pwr dp ' + 7- m 

A further integration of (|30|) followed by a multiplication with p yields 

= --+P f 4* f 2p"|y(p") + 2pV(p")| 2 ^"^' + l + C 2 p. (31) 



»■ ' ' " 7o p' 2 ./o 

As in [5] we find that C2 = k^(0) which is usually, but not always, put to 0. 
This gives lemma |3j □ 

Using this lemma it is not difficult to prove that y(p) and thus m n specifies 
/i. In the next section we will see that the field equations can be solved with 
the condition g(r, 0) = /i(r). 



Theorem 4. Given a set of multipole moments m n such that (|22p has posi- 
tive radius of convergence, the corresponding fi, is uniquely specified and real- 
analytic in a neighbourhood of the origin. 

Proof. The proof is essentially the same as in [2]. From (f2"3"|) we find that 

r(p)=p/(l + p S (p)), (32) 

where 

sOO = f 4s f V Mp") + VWOI 2 dP 7 + 1 + K '(o)- (33) 

Jo P Jo 

is real-analytic in a neighbourhood of p = 0. From r(0) = 0, r'(0) ^ 0, the 
inverse function theorem for analytic functions gives that p = p(r) is also a 
real-analytic function, with p(0) = 0, in some neighbourhood of r = 0. The 
equation (|24|) can be written 



hi?) = = ^l + p(r)Q(p{r))y{p(r)). (34) 

Thus /l (r) is a complex valued real-analytic function in a neighbourhood of the 
origin. □ 
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6.1 Solving the field equations 

We can now turn our attention to the problem of extending /l to a solution / 
of the field equations. 

Lemma 5. Let be analytic in a neighbourhood of the origin. Then there 

is a formal power series solution C) to (|14|) such that g(£,0) = 

Proof. In a neighbourhood of the origin we can formally expand the non-linear 
term of f| 14[) . Thus, define c n m via 



oo oo 



E E ^cr - 2 ^ + ^ + 2 ^\ (35) 



Due to the symmetry in £ and £, we have Cn jm = c m ,n- 

Let P„ denote the set of polynomials in ak,i and a^ i with non-negative 
coefficients, where k + I < n. We show that c„ i?Ti = c nim (a. t .,a. t .) S V n +m- The 
£ n £ m -coefficient of (5 + 2^g^)(g + 2£<jtf) is easily found to be 

n m 

E E «»n-fc,lOm-l,fc(l + 2(n - fc))(l + 2(m - /)). (36) 

fc=0 2=0 

Thus it is a polynomial in dk,i with non-negative coefficients, where k + I < 
n + m. The £ n £ m -coefficient of (gg) p , where p € N is easily seen to be in V n+m . 
Therefore, the ^"^-coefficient of 

1 00 

T^rX>«»r (37) 

is also in Vn+m- Consider the set of power series in £ and £, with the £ n £ m - 
coefficient in V n +m- This set is closed under products. Thus c„. m G V n + m - 
The equation (|14|) can be rewritten as the system 

fe k + n 

(k + l)(n + |)a n+ i ife+ i = E (! - m H 2 — ) a ' l + fe + 2 -"^" 1 - c ». fc ( 38 ) 

m— 

where < k < n. 

This then can be reformulated as 

(2A; - l)(n - fc + 1) c n - k ,k-2 - Cn-k-i,k-i , on x 

a -^ fc " k{2n-2k+l) "-W- 1 + 2 fc(2n-2A ; + l) (39) 

and simplified to 

(2fc - l)!!(2n - 2fc - 1)!! / n!a„, 2(fc - l)!(n - k)\cn- k -i,k-i 



a n -k,k 



k\{n-k)\ \(2n-l)\\ (2fc-l)!!(2n-2jfe + l)!! 



fe— 1 



(z - l)!(n - i - l)!(n - 2i) 



^ (2i+l)!!(2n-2« + l)! 



"Cn— 2 — 1,2 — 1 



(40) 



where 1 < k < [f J ■ The data <?(£, 0) gives a n ,o and together with (J3S]), the 
recursion (|40[) gives a formal power series solution £). □ 
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Now we want to establish convergence of the series produced by the previous 
lemma. For this we will use the technique of majorizing series. 

Definition 6. A series X)^°=o ^ n b n > is said to majorize a series 

S^°=o a " w tth a n > if b n > a n for all n. 

Lemma 7. The formal solution from lemma^ converges in a neighbourhood of 
the origin. 



Proof. Now let M,R > be such that —M — majorizes J2n=o \ a n,o\£, n - This 

can always be done. For instance choose R such that the series is bounded on 
< £ < 2i?, and let M be the maximal value of the series on that interval. 

The triangle inequality on the recursive step (|40p . together with the non- 
negativity of the coefficients of c, implies 

. (2k - l)!!(2n -2k- 1)!! ( n\\a n . \ , 2{k - l)!(n - fc)!|cn-fc-i,fc-i| 

|On-fc,fc| S 

fc-1 



k\(n-k)\ \(2n-l)\\ (2k - l)!!(2n - 2k + 1)!! 



< 



(i-l)!(n-i-l)!(n-2i) , 
^ (2i + l)!!(2n-2i + l)l! |c "- 1 - 1 ^ 11 

(2k- l)!!(2n-2fc- 1)!! / n!|a„, | 



fc!(n-ft)! Vl 2 "- 1 )- 
2(fe - l)!(n - fc)!cn-fc-i,k-i(k,.|, |a.,.|) 



(2fc - l)!!(2n - 2fc + 1)!! 



fe-i 



„^ (t-l)!(rc-i-l)!(n-2z) 

(2i + l)!!(2n-2i + l)!! c — M-i(KI, kJ)^ 



(41) 



Now we want to construct a series majorizing X)^o Em=o \ a n,m\£, n ( m ■ Let 
bn,fc > satisfy 

_ (2k- l)!!(2n-2fe- 1)!! / n! 
"- fc ' fc " fc!(n-A;)! ^ (2n - 1)!! "'° 

2(fc- l)!(n-fc)!c n _ fc _ 1 ^ 1 (b.,.,b.J 

(2A-l)!!(2n-2Jfe + l)!! ( 42 ) 

i0 ^ (z-l)!(n-z-l)!(n-2») \ 
+ 2 Z. (2i + l)!!(2n-2i + l)!! c — 

and form 

OC 71 — 1 

n— Tci— 

By comparing with lemma [51 we see that recursion (|42[) gives the coefhcients of 
the analytic solutions of the equation 

% = 2iw) + r-icc^ ■ (44) 
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This equation can be simplified to a linear equation by a non-linear transforma- 
tion. To find this transformation we substitute 



f(V^«(e,o) j (45) 



in (I44p and force the a^a^-term to vanish. The resulting differential equation 
for f can then be solved to obtain the transformation 

sinfi arcsin(3\/£CW£, C))) , > 

hfo = — ^ (46) 



Substituting this into (f44|) . we get the linear equation 

"« - m^tr m 

If we let a(£, 0) = 0) = , M c , we get the solution a(£, C) - 
thus 



sin(| arcsin( 3A/ ^ )) 

fc(f,0 = J^c ™ 

Note that this is real-analytic in £ and C, with non-negative coefficients, and 
have the symmetry h(£, £) = Furthermore it is bounded in the region 

u — ' "s — 1+3MB. ■ 

The choice of R,M gives |a nj o| < &«,o- We will now prove that |a„_fc^| < 
b n -k,k for all < k < [§J and all n £ N by induction over n. The base step 
(n = 0, 1) is trivial. The induction hypothesis is that it holds for n — 2. Together 
with the inequalities |a n o| < ^n,o an d (14ip we have 



(2fc-l)!!(2n-2fc-l)!! / n! 
K - fc ' fc| - k\{n-k)\ U2n-1)!! n '° 

2(fc-l)!(n-A:)!c n - fc _i, fc -i(b.,.,b.,.) 

(2fc-l)!!(2n-2fc + l)!! (4Q) 

^ (i-l)!(n-i-l)!(rc-2i) A 
+ 2 Z. (2i + l)!!(2n-2i + l)!! ' " 1 ; J 

bn—k,kj 

where the last equality is Hence |a n ,fc| < &n,fc for arbitrary n, fc £ N with 

k < n. Thus the series expansion for 5 converges whenever h does. □ 

Note that by lemma [5] we can compute a number of terms of the power series 
g. The modulus of the remaining terms can then be estimated by the methods 
in lemma [7] However for the proof of existence of solutions we only need the 
convergence. 

Now we can map <?(£, £) back to f(w) and change back to the variables (z, p) 
via w = Z+ 2 P ■ The potential will still be analytic and even in p. Therefore we 
have analyticity in z and p 2 = x 2 + y 2 . Thus, together with theorem 0] we have 



10 



Theorem 8. For every set of multipole moments m n such that the power series 
(I22[) has positive radius of convergence, there is a unique solution f to the field 
equations, in the class of analytic functions of {x,y,z). 

Now we have the potential /, we easily find A and lo. It can be shown that 
the metric function j3 can always be obtained by integrating (|20p . The constant 
of integration can always be chosen such that &(£,C) — KC>£) an d that ^t'^yi 
is analytic in £ and £. With this choice the corresponding (3 will be analytic in 
(x,y, z), and tend to zero fast enough for V to be asymptotically flat. Hence, the 
solution has all properties demanded in the theorems. Therefore the multipole 
moments can be freely specified as long as the power series (|2"2"j) has positive 
radius of convergence. 

The metric function W can also be obtained by an integration. If we don't 
demand mo to be real W will not asymptotically tend to zero. In this case we 
get a non-zero NUT-parameter. 

7 Conclusions and discussion 

We can conclude that if the multipole moments are chosen such that the power 
series ([22"|) has positive radius of convergence, the corresponding potential func- 
tion / exist and is analytic in the rescaled spatial variables. Furthermore the 
corresponding metric functions for the spacetime can be found. Asymptotic 
flatness of the manifold of time-like Killing trajectories is guaranteed, but the 
spacetime will in general not be "asymptotically Minkowski". 

The methods presented here can also be used to construct axisymmetric 
spacetimes with specified multipole moments. The leading order part of the 
potential can be found by theorem^ However this involves a function implicitly 
defined, which can make explicit computations difficult. The main difficulty is 
to solve the recursion in lemma [51 The recursion can however be used to get 
approximations to arbitrary order. The methods in lemma [7] can then be used 
to get estimates of the solution or the rest terms. 

For the future it is natural to try to extend the results to the general station- 
ary case, without the axisymmetry condition. The field equations will then be 
much more complicated, therefore this task will be difficult. Another possibility 
is to extend the results of this paper to the electro-vacuum case. 
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